Unusual behaviour of the ferroelectric polarization in PbTiOs/SrTiOs superlattices 
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Artificial PbTiOs/SrTiOs superlattices were constructed using off-axis RF magnetron sputtering. 
X-ray diffraction and piezoelectric atomic force microscopy were used to study the evolution of 
the ferroelectric polarization as the ratio of PbTiOa to SrTiOa was changed. For PbTiOa layer 
thicknesses larger than the 3-unit cells SrTiOa thickness used in the structure, the polarization is 
found to be reduced as the PbTiOs thickness is decreased. This observation confirms the primary 
role of the depolarization field in the polarization reduction in thin films. For the samples with 
ratios of PbTiOa to SrTiOa of less than one a surprising recovery of ferroelectricity that cannot be 
explained by electrostatic considerations was observed. 



The construction of artificial ferroelectric oxide super- 
lattices with fine periodicity presents exciting possibili- 
ties for the development of new materials with extraor- 
dinary properties and furthermore is an ideal probe for 
understanding the fundamental physics of ferroelectric 
materials. 

The most studied system at present is BaTiOa/SrTiOs 

Other combinations that 
have been studied include KNbOa/KTaOa H 
El E3> PbTiOa/SrTiOa [H, PbTiOs/BaTiOa „ 
PbTiOs/PbZrOa'llEl and most recently high qu ality 
tricolour superlattices of SrTiOs/BaTiOg/CaTiOslH. 
In BaTiOs/SrTiOa, first principles studies suggest 
that both the SrTiOs and BaTiOs layers are polar- 
ized such that the polarization is approximately uniform 
throughout the superlattice. The driving force behind 
this is the large electrostatic energy penalty for a build- 
up of charge at the interface caused by discontinuous 
polarization in the normal direction. The electrostatic 
model proposed by Neaton and Rabe |5i| to explain their 
first principles results for BaTiOs/SrTiOa superlattices 
is very similar to the electrostatic model applied to cal- 
culate the effect of the depolarization field in ultra-thin 
ferroelectric films with realistic electrodes 0, 0, Eoj . 
Experimentally it was recently shown that the reduced 
polarization observed in monodomain thin PbTiOa can 
be explained by the presence of a depolarization field re- 
sulting from imperfect screening of the polarization (ij . 
Recent work also suggests that, under certain conditions, 
the electrostatic energy due to depolarization fields will 
drive the system to form domains as observed by Fong et 
al. '2^ and Nagarajan et al. In this letter we use 

PbTiOa/SrTiOa superlattices to probe the effect of a re- 
duced ferroelectric thickness in a dielectric environment. 
Our data show that the behaviour observed in PbTiOa 
thin films is reproduced for PbTiOa layers thicker than 
three unit cells. However, for thinner ferroelectric layers 
a surprising recovery of ferroelectricity that cannot be 



explained by electrostatic considerations was observed. 

The superlattices of PbTiOa/SrTiOs were prepared 
on conducting 0.5% Nb doped (001) SrTiOs substrates 
using off-axis RF magnetron sputtering with conditions 
similar to those used for growing high quality epitax- 
ial c-axis PbTiOa thin films For all the samples 
discussed in this paper, the SrTiOa thickness was fixed 
at three unit cells (about 12 A). At room temperature 
the in-plane lattice parameters of tetragonal ferroelec- 
tric PbTiOa (a=3.904 A,c=4.152 A) and cubic dielectric 
SrTiOs (3.905 A) are an excellent matchjl^. It is thus 
expected that PbTiOa will grow coherently on SrTiOa 
substrates, and that the strain interactions will be dom- 
inated by the constraint imposed by the substrate. The 
growth temperature for the superlattices was 460° C. In- 
vestigation by Transmission Electron Microscopy (TEM) 
revealed excellent quality in superlattices with layers of 
SrTiOa thinner than 5 unit cells, though beyond this 
thickness the quality of the SrTiOs layers deteriorated 
with thickness, presumably because of the low tempera- 
ture. On the other hand, samples processed with higher 
temperatures were of lower quality, probably because 
of lead losses from the PbTiOa. The low temperature 
growth used thus seemed to be optimal for making su- 
perlattices in which the SrTiOa layers are very thin, but 
without limitation on the thickness of PbTiOa. 

In the principal series of interest we grew superlattices 
consisting of 20 PbTiOs/SrTiOs bilayers in which the 
SrTiOa layer thickness was maintained at 3 unit cells 
while the PbTiOs layer thickness n was varied from 54 
unit cells down to just one unit cell (denoted n/3). The 
first layer deposited was PbTiOs. The layer thicknesses 
were determined by calibrating x-ray diffraction patterns 
with deposition time. 

Cross-sectional TEM investigations were performed on 
several samples and reveal the coherent growth and ar- 
tificial layering of the samples. Fig. ^ shows a sum- 
mary of the results obtained on a 3/3 sample. The bright 
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field image, Fign(a), siiows tiie layering throughout the 
sample. The perfect crystalline structure and coherent 
growth are demonstrated by the high resolution TEM 
(HRTEM) image Fig(T](b), while the periodicity of the 
superlattice is demonstrated by the superlattice reflec- 
tions in the diffraction image Fign(c) (arrows). 




FIG. 1: Cross sectional TEM images of a 20 bilayer 
PbTiOz/SrTiOi 3/3 sample, (a) Bright field image clearly 
shows the intended layering of the structure, (b) HRTEM 
shows the perfect crystalline structure of the material, (c) 
Diffraction image demonstrating superlattice periodicity. 

Further structural characterization were performed us- 
ing standard 9 — 26 x-ray diffraction. Fig. [21 shows the 
— 26 scan for a superlattice in which the PbTiOa layers 
are 9 unit cells thick and the SrTiOs layers are 3 unit cells 
thick (9/3). The periodicity of the superlattice is there- 
fore 12 perovskite unit cells and 12 reflections from 29 — 
to the angle corresponding to the 001 peak of the average 
perovskite unit cell lattice parameter (at 26 « 22°) are 
expected, most of which are observable in the scan. In 
between the main superlattice peaks, the presence of 18 
finite size effect peaks, clearly visible in the inset of Fig. 
El is due to the finite total thickness of the sample which 
is 20 times the superlattice periodicity. 




FIG. 2: 6 — 26 x-ray diffractogram for a 20 bilayer 
PbTiOs/SrTiOi 9/3 superlattice. 

Because of the large strain-polarization coupling in 
PbTiOa [i^l , a change in polarization results in a change 



in material tetragonality |2l| . We take advantage of this 
to follow the evolution of the polarization in the superlat- 
tice by following the evolution of the average c axis lattice 
parameter, c, as the PbTiOa layer thickness is varied. If 
the wavelength of the superlattice is nc then the n*'' peak 
Yo. a, 9 — 29 scan will always correspond to c irrespective 
of the value of n allowing the average c axis lattice pa- 
rameter of the superlattice to be easily determined. In 
practice, this peak is easily identifiable due to its high in- 
tensity and proximity to the substrate peak. Intuitively 
one expects, as the thickness of the PbTiOa layers rel- 
ative to the SrTiOa layers is reduced, a decrease of the 
ferroelectric polarization which should result in a con- 
comitant decrease of the average lattice parameter. The 
measured average c axis lattice parameters as a function 
of the thickness of the PbTiOa layer thickness are shown 
in Fig. O For comparison, we also show the average c 
axis lattice parameters obtained by fixing c of SrTiOs 
at its paraelectric cubic value 3.905 A and taking c of 
PbTiOa in two limiting cases: first, at the value 4.022 A 
corresponding to a hypothetical paraelectric tetragonal 
structure coherent with the substrate (solid line)|2l| and 
then at the fully polarized bulk value 4.152 A (dashed 
line). As can be seen in Fig. (3 superlattices with thick 
PbTiOa layers have "large" average lattice parameters 
clearly suggesting a ferroelectric polarization. On reduc- 
tion of the layer thickness the average lattice parameter 
decreases and approaches the solid line. However, sur- 
prisingly, after reaching this line superlattices with very 
small PbTiOs layer thicknesses display larger average lat- 
tice parameters which indicate a recovery of ferroelectric- 
ity. 

This behavior was confirmed using atomic force mi- 
croscopy (AFM) which allows the ferroelectric domain 
structure to be modified and detected on a local scale jl^. 
Applying a voltage between the metallic tip of the AFM 
and the metallic substrate, stripes were "written" (poled) 
using alternatively positive and negative voltages. Piezo- 
electric atomic force microscopy (PFM) was then used 
to detect the domain structure. PFM images are shown 
in the insets of Fig. O for different superlattices, the 
contrast revealing domains with up and down polariza- 
tion. As can be seen, the 1/3, 2/3 and 13/3 samples 
reveal a clear domain structure and are indeed ferroelec- 
tric whereas no significant contrast could be obtained in 
the 3/3 superlattice, confirming the behavior suggested 
by the x-ray analysis. The written domains for all sam- 
ples in which domains could be written were confirmed 
to be stable for a number of days. All domains writ- 
ten into the down direction have the same piezoelectric 
response as the existing background, implying that be- 
fore writing the entire sample is uniformly poled in the 
down direction, demonstrating that none of the samples 
formed a polydomain state. This is a particularly im- 
portant observation for the 1/3 and 2/3 samples since 
a multidomain configuration could possibly explain the 
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observed increase in average c axis lattice parameter at 
very small PbTiOa layer thicknesses. 
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FIG. 3: Average c-axis lattice parameter plotted against the 
number of unit cells of lead titanate per bilayer showing the 
suppression and recovery of ferroelectricity. Complementary 
PFM images are shown as insets. 

To understand the observed behaviour, a simple elec- 
trostatic model similar to the one proposed by Junquera 
and Ghosez has been developed. The total energy 
per unit cell area E of an np/ris super lattice is written 
as 

= Up Up{P^) + 71, C/,(P0) + Peiec(P;, P."), 

(1) 

where Up and Us are the total energies per 5-atom unit 
cell of bulk PbTiOs and SrTiOa in zero field as a func- 
tion of their polarization Pp and P° (assumed to be ho- 
mogeneous in each layer) and E^iec is the macroscopic 
electrostatic energy resulting from the presence of non- 
vanishing electric fields in the layers when P^ and P° 
differ. 

The electrostatic energy of a given layer, of thickness 
Ip or Zs, in the presence of a finite electric field to 
leading order in the field, is Eeiec = —^^ ■ In the 
superlattice, the electric fields £p and Es are determined 
by Pp and P° through the condition of continuity of the 
normal component of the electric displacement field at 
the interfaces; 



Pp + Eo^p — Ps + £oEs 



(2) 



For a system under short-circuit boundary conditions, 
the potential drop along the structure must vanish so 
that 



IF — -IF 



(3) 



Combining the last two conditions and summing the elec- 
trostatic energies of the PbTiOa and SrTiOa layers we 
obtain 



i?e/ec(P°,P°) 



eo{lp + ls) ' 



The total energies Up{Pp) and Us{Ps) have been ob- 
tained from density functional theory (DFT) calcula- 
tions on bulk compounds using the ABINIT package |27| . 



The calculations were performed within the local den- 
sity approximation (LDA) using extended norm con- 
serving pseudopotentials [29| with the Pb (5d,6s,6p), Sr 
(4s,4p,5s), Ti (3s, 3p, 3d, 4s) and O (2s, 2p) treated as va- 
lence states. Convergence was reached for a 1225 eV (45 
Ha) cutoff and a 6 x 6 x 6 mesh of special fc-points. We 
obtain computed lattice constants for cubic paraelectric 
SrTiOa (a = 3.846 A) and for tetragonal ferroelectric 
PbTiOa (a = 3.864 A, c = 3.975 A), with a polariza- 
tion of 0.69 C/m^. Both materials lattice parameters 
are underestimated relative to the experimental values 
as is typical for the LDA. For each compound, U{P^) 
and c(P^) were obtained [s^ following the formalism of 
Ref. |28| by relaxing the atomic positions and the lattice 
parameter c at fixed polarization P^z in the space group 
PAmm, constraining the in-plane lattice parameter a to 
3.846 A. For bulk PbTiOs constrained in plane to 3.846 
A the c axis lattice parameter was found to be 4.009 A 
with a polarization of 0.73 C/m^. For any choice of Up 
and ris , minimization of Eq. ^ gives equilibrium values 
for Pp and P", and thus also for Cp and Cg- 

To see whether the model correctly describes the be- 
haviour as the thickness of the PbTiOa layers decreases 
to the atomic scale, we performed full DFT-LDA calcula- 
tions of the structure and polarization of PbTiOa/SrTiOa 
superlattices for = 3 and rip = 1,... 7 with the Vienna 
ab initio Simulations Package (vASP) 1321 . using projec- 
tor augmented wave (PAW) potentials with the 
same valence configurations as in the ABINIT calculation. 
Convergence was reached for a 600 eV (22 Ha) cutoff and 
a 6 X 6 X 2 mesh of special /c-points. The computed lattice 
constants are for SrTiOs, a — 3.86 A, and for tetragonal 
ferroelectric PbTiOs, a = 3.86 A and c = 4.047 A, with 
a polarization of 0.75 C/m^ [s^. For the superlattices, 
the atomic positions and lattice parameter c were fully 
relaxed in the space group PAmm, constraining the in- 
plane lattice parameter a to 3.86 A. Polarizations were 
calculated using the modern theory of polarization [ssj 
as implemented in VASP. 

Fig. 0] (a) shows the evolution of the polarization as a 
function of ^ for = 3. According both to the model 
and to the first principles local polarizations (not shown), 
the difference between the polarizations in the two layers 
is quite small, highlighting the large electrostatic energy 
cost of having different polarizations in the layers. As the 
ratio Up/us increases, the polarization of the superlattice 
asymptotically approaches the constrained bulk PbTiOs 
value, though rather slowly due to the large energy cost 
of maintaining a high polarization in SrTiOa . The figure 
inset shows the corresponding increase in the tetragonal- 
ity (c/a) of the two layers, with the high polarization- 
strain coupling in the SrTiOa layer (higher even than for 
the PbTiOa layer l^^l) being evident. Both the model 
and the first principles calculations show a monotonic 
decrease of the polarization as the PbTiOs volume frac- 
tion is reduced, due to the increase in the relative energy 
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cost of the polarization in the SrTiOs layers. While the 
polarization vs thickness curve for the model at the low- 
est thicknesses is shifted to lower polarizations relative to 
the first principles results, the model works overall very 
well, considering the simplifying assumptions and lack of 
any adjustable parameters. 




a 4 « e 10 12 14 16 
iip/n. 



FIG. 4: (a) Polarization in each layer obtained from the elec- 
trostatic model (dotted and solid lines) and the average po- 
larization from first principles calculations (full circles). In- 
set shows tetragonality in each material calculated from both 
methods, (b) Comparison of experiment and both theoretical 
approaches. 

In Fig. ^ (b) we compare results from the first prin- 
ciples calculations (open circles), the electrostatic model 
(solid line) and experiment (solid squares), by plotting the 
fractional change in the superlattice tetragonality f — 1 
relative to the tetragonality of bulk PbTiOa with the in 
plane lattice parameter constrained to the SrTiOs sub- 
strate 1^ . Good agreement between both theoretical ap- 
proaches and experiment is seen for samples that are pre- 
dominantly PbTiOa. It should be noted that both the- 
oretical calculations are at zero temperature, whilst the 
experiments are conducted at room temperature. Specif- 
ically this means that samples predicted from first prin- 
ciples to be ferroelectric with a small polarization at zero 
temperature might be expected to be paraelectric in our 
room temperature experiment, as is observed in the case 
of the 3/3 sample. The fact that unexpected recovery 
of the ferroelectric polarization in the experimental 1/3 
and 2/3 superlattices is observed in neither the electro- 
static model, nor the first principles calculations, sug- 
gests that it is related to aspects not accounted for in 
our theoretical approaches, for example, the precise na- 



ture of the substrate-superlattice interface, some degree 
of intermixing at the superlattice interfaces, or the pos- 
sible formation of a new entropically stabilized PbTiOs 
phase similar to that formed under negative hydrostatic 
pressure in the first principles studies of Tinte et al [3^ . 
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